Evaluation of quadrature schemes for the discrete ordinates method

After reviewing the presently available quadrature schemes for the discrete ordinates method, the accuracy of different schemes is analyzed and evaluated. It is shown from a comprehensive error analysis that the moment conditions have to satisfied not only for the principal coordinates directions, as it is mostly carried out, but for any arbitrary test direction. Among the schemes with approximately 50 discrete ordinates the DCT-020-2468 quadrature was found to give the best accuracy. The highest accuracy among all schemes is achieved by the LC-11 quadrature which requires 96 discrete ordinates. This scheme is rarely used up to date and deserves more attention for high accuracy predictions.     

Plastic multipliers as driving variables of numerical simulation in elastoplasticity

Statics of elasto-plastic media is stated in terms of eigenstrains acting upon a background elastic problem with fixed (initial) stiffness. In order to minimize the number of unknowns and to provide computationally cheap algorithms plastic multipliers are used as main driving variables. Fixed-point type iterations are suggested for computing plastic multipliers within the load increment. Numerical experiments results are shown to be in good correlation with other software results. A range of problems is outlined where the usage of the proposed algorithms can be advantageous.     

Large-scale computer simulation of fully developed turbulent channel flow with heat transfer

Recently, with the advent of supercomputers, there has been considerable interest in the use of direct numerical simulation to obtain information about turbulent shear flow at low Reynolds number. This paper presents a pseudospectral technique to solve the full three-dimensional time-dependent Navier-Stokes and advection-diffusion equations without the use of subgrid-scale modelling. The technique has not been previously used for fully developed turbulent channel flow simulation and is based on methods applied in other contexts. The emphasis of this paper is to provide a reasonably detailed account of how the simulation is done rather than to present new calculations of turbulence. The details of an algorithm for turbulent channel flow simulation and the grid and time step sizes needed to integrate through transient behaviour to steady state turbulence have not been published before and are presented here. Results from a Cray-2 simulation of fully developed turbulent flow in a channel with heat transfer are presented along with a critical comparison between experiment and computation. The first- and second-order moments agree well with experimental measurements; the agreement is poor for higher-order moments such as the skewness and flatness near the walls of the channel. Detailed information given about the effects of spatial grid resolution on a computed results is important for estimating the size of the computation required to study various aspects of a turbulent flow.     

FULL DISCRETE NONLINEAR GALERKIN METHOD FOR THE NAVIER-STOKES EQUATIONS

This paper deals with the inertial manifold and the approximate inertial manifold concepts of the Navier-Stokes equations with nonhomogeneous boundary conditions and inertial algorithm. Furtheremore, we provide the error estimates of the approximate solutions of the Navier-Stokes Equations.     

Standing waves, clustering, and phase waves in 1D simulations of kinetic relaxation oscillations in NO+NH3 on Pt(1 0 0) coupled by diffusion

The Lombardo–Imbihl–Fink (LFI) ODE model of the NO+NH₃ reaction on a Pt(1 0 0) surface shows stable relaxation oscillations with very sharp transitions for temperatures T between 404 and 433 K. Here we study numerically the effect of linear diffusive coupling of these oscillators in one spatial dimension. Depending on the parameters and initial conditions we find a rich variety of spatio-temporal patterns which we group into four main regimes: bulk oscillations (BOs), standing waves (SW), phase clusters (PC), and phase waves (PW). Two key ingredients for SW and PC are identified, namely the relaxation type of the ODE oscillations and a nonlocal (and nonglobal) coupling due to relatively fast diffusion of the kinetically slaved variables NH₃ and H. In particular, the latter replaces the global coupling through the gas phase used to obtain SW and PC in models of related surface reactions. The PW exist only under the assumption of (relatively) slow diffusion of NH₃ and H.     

Numerical analysis of roughness effect on microtube heat transfer

Roughness effect on the heat transfer and pressure loss performances of microscale tubes and channels are investigated through a finite element CFD code. Surface roughness is explicitly modelled through a set of random generated peaks along the ideal smooth surface. Different peak shapes and distributions are considered; geometrical parameters are representative of tubes in the diameter range from 50 to 150 μm. The use of a fine enough mesh allows the direct computation of tube performances under the assumption of incompressible fully developed flow. As a result, a significant increase in Poiseuille number is detected for all of the configurations considered, while the effect of roughness on heat transfer rate is smaller and highly dependent on the tube shape.     

分离变量法与等腰直角三角形波导的边值问题

根据叠加原理，讨论等腰直角三角形波导边值问题的IE波解和TM波解，纠正有关文献中的错误解法，指出分离变量法的适用条件。     

Full-Wave Analyses for a Multi-Layered Microstrip Patch Antenna

This paper describes numerical modelling of a dual band multi-layered microstrip patch antenna operating at 35GHz on an in-house semi-insulating GaAs substrate. The simulated and measured resonant frequencies at both lower and upper resonant frequencies will be compared to check the accuracy of the different numerical modelling techniques.     

Comparison between Eulerian and Lagrangian semi-analytical models to simulate the pollutant dispersion in the PBL

In this work we display a numerical comparison, under statistical and computational point of view, between semi-analytical Eulerian and Lagrangian dispersion models to simulate the ground-level concentration values of a passive pollutant released from a low height source. The Eulerian approach is based on the solution of the advection–diffusion equation by the Laplace transform technique. The Lagrangian approach is based on solution of the Langevin equation through the Picard’s Iterative Method. Turbulence inputs are calculated according to a parameterization capable of generating continuous values in all stability conditions and in all heights of the Planetary Boundary Layer (PBL). Numerical simulations and comparisons show a good agreement between predicted and observed concentrations values. The comparison reveals the main advantages and disadvantages between the models.